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ABSTRACT 

I n this paper £he crack problem for a relatively thin layer bonded 
to a very thick substrate under' "thermal shock" conditions is considered. 
The effect of surface cooling rate is studied by assuming the temperature 
boundary condition to be a ramp function. Among the crack geometries 
considered are the edge crack in the coating layer, the broken layer, the 
edge crack going through the interface, the undercoat crack in the sub- 
strate and the embedded crack crossing the interface. The primary calcu- 
lated quantity is the stress intensity factor at various singular points 
and the main variables are the relative sizes and locations of cracks, 
the time, and the duration of the cooling ramp. The problem is solved 
and rather extensive results are given for two material pairs, namely 
a stainless steel layer welded on a ferritic medium and a ceramic coat- 
ing on a steel substrate. 

1. Introduction 


A common failure mode in many structural components that consist of 
a relatively thin coating and a substrate is the cracking of the coating, 
the substrate, or both under transient thermal stresses. Cladded pres- 
sure vessels under thermal shock, and certain microelectronic devices and 
ceramic-coated metal parts under rapidly changing thermal environments 
may be mentioned as some typical examples for such components. Very 
often in transient thermal stress analysis it is assumed that the rele- 
vant thermal boundary condition is a step change in temperature. In the 
fracture analysis even though this assumption would lead to conservative 
bounds for the stress intensity factors, in most cases it is not a real- 
istic representation of the actual boundary conditions. 

The broad aim of this study is to provide the solution needed for 
assessing the crack propagation and arrest process in a coated medium 
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containing a certain initial flaw and subjected to transient thermal 
loading. In a simplified fracture analysis this requires the determina- 
tion of the stress intensity factors as functions of time and dimensions 
concerning the size and location of the crack. In order to study the 
influence of the thermal boundary conditions on the fracture process, in 
this study the rise time of the surface temperature is assumed to be an 
additional variable in formulating the crack problem. 

The crack problem is solved under the following simplifying assump- 
tions: (i) the medium consists of two dissimilar linear isotropic homo- 

geneous materials bonded along an ideal plane interface; (ii) the crack 
lies in xz plane, the constraints in z direction are such that the prob- 
lem may be treated as either a plane strain or a generalized plane stress 
problem in xy plane, and the thermal boundary conditions are independent 
of y; (iii) the thickness of the substrate is very large compared to that 
of the coating, hence for small times the main features of the problem 
may be recovered by approximating the medium by a layer bonded to an 
elastic half space; (iv) all thermoelastic coupling effects and the tem- 
perature dependence of thermoelastic coefficients are negligible; and 
(v) the transient thermal stress problem is quasi static, that is the 
inertia effects are negligible. The previous studies on dynamic thermo- 
elasticity seem to bear out the validity of the last assumption (see, 
for example [1] and [2]). 

The thermal shock problem for hollow cylinders containing an axial 
or a circumferential crack was considered in [3]-[5], where it was 
assumed that the thermal boundary condition is a step function in tem- 
perature. In 13] the axisymmetric circumferential crack problem was 
analyzed for a homogeneous hollow cylinder. The same problem was studied 
in [4] by assuming that the cylinder is homogeneous in elastic but non- 
homogeneous in thermal properties, simulating the cladded pressure ves- 
sels. In [5] the crack is assumed to be axial, that is, in the rz plane 
and the coating is approximated by a "membrane" in the elasticity solu- 
tion and is assumed to be a homogeneous continuum in solving the diffu- 
sion problem. 
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2. Temperature Distribution 

Consider the nonhomogeneous thermoelastic medium shown in Fig. la. 
Assume that initially the medium is at a constant temperature To, and 
starting at t=0 the temperature of the surface x=0 is (suddenly or 
gradually) changed to T Q and is held constant at T Q for t>t 0 (Fig. lb). 
Defining 

Tj(x»t) - Ta> = e|(x,t) , i=l ,2 , (1) 

the differential equations and the boundary conditions for the heat con- 

(*) 

duction problem may be expressed as' ' 


a z e i 

3X 2 


D i 


36 . 

HT 


(1=1,2) 


( 2 ) 


e^x.O) = 0, (i=l ,2) , 

39,(h,t) 36 ? (h,t) 

e^hjt) = 0 2 (h,t), — = k 2 ' — , 


(3) 

(4) 


02 ( oo » t ) = 0 


(5) 


e-j(0,t) = e Q H(t) , e Q = T 0 ~T 00 , or (6) 

© 1 (0, t) = ^ [tH(t)-(t-t Q )H(t-t 0 )] , e Q = T 0 -T ro . (6a) 

u o 

The solution of (2) subject to (3)— (6) may be obtained in a straightfor- 
ward manner by using Laplace transforms (see, for example, I6],[7j). 
Defining the dimensionless quantities 


x 1 = x/h, 6 = /D^7Dj , n = 6k 1 Vk 2 ' , t = tD^h 2 , (7) 

for the boundary condition (6) the temperature distributions in the two 
materials are found to be 


e^xV) 



f -TS 2 . | 

e sinx's d 
s(cos z s+n z sin z s) 

0 


(*) See the List of Symbols for notation 


0<x'<l , 


( 8 ) 
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9o(x',t) 2 r°° 2 

-=— = 1 -- {e s [cos s sine(x ' -1 ) s 


+ nsin s cosb(x'- 1 )s]/[s(cos 2 s+ri 2 sin 2 * s)]}ds, l<x'«» . (9) 

Similarly, if the temperature boundary condition at x=0 is given 
by the ramp function shown by (6a), we obtain 

e l(x ,t) 2 n f (l-e“ Ts2 )sinx's ds n , fW'r Mf)> 

§ ’ ° <X ' 5 ° T T o ’ (10) 

n n n * x 


9 1 (x,,t) = 1 _2n 


7TT 


' %Ssw;l&r' s ds - ° <x,<l - t>t o <"> 


^2^ » T ) T 2 

e 


T 7TT 

0 O 


-TS 2 


{(1-e )[cos s sine(x'-l)s 


+ nsin s cos8(x , -l)]/[s 3 (cos 2 s+n 2 sin 2 s)]}ds, x'>l, 0 <t<t , 


( 12 ) 


©2( x 2 

8 7TT 

O O 


r c 2 ^ s z 

{e" TS (e 0 -1 )[cos s sim3(x'-l )s 


+ nsin s cosb(x'- 1 )s]/[s 3 (cos 2 s+ri 2 sin 2 s)]}ds, x’>l, t>t 


where 


T o = W h2 


(13) 

(14) 


3, Thermal Stresses in the Uncracked Medium 

In the uncracked medium subjected to thermal initial and boundary 
conditions (3)-(6) x and t are the only independent variables, the medium 

is unconstrained in x direction and is fully constrained in y and z 
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directions. For the thermal stresses we thus have 

°ixx ~ ° 5 e ixy ” e ixz ~ £ iyz " e iyy ” G izz “ a iyy _CT izz’ ^ 4 5 * * * ) 

From (15) it follows that 


T 

°iyy 


(x,t) 


T F-|0t-l ®1 t) 

a lz2 (x,t) - - T~v 5 0<x<h, t>0 , 


(16) 


a 2yy^ x ’^ ^zz (x »t) 


: 2 a 2 0 2 (x 
“1 ^ 


, x>h, t>0 . 


(17) 


By defining 


_ E l a l 9 o , _ x tD l 
a o 1-v-j * x ‘ h » T ~W 


(18) 


the stresses may also be expressed in the following normalized form: 



9-] ( X »T ) 

e 

0 


0 < x ' <1 , t >0 , 


p 2yy (x '» T) 

a 

o 


E 2 a 2 (l-v 1 ) 6 2 (x',t) 
E 1 a 1 (l-v 2 ) 6 o 


X 1 > 1 , T>0 . 


(19) 


(20) 


4. The Crack Problem 

Consider now the crack problem shown in Fig. la. Since the problem 

is linear, one may consider the solutions due to thermal and mechanical 

loads separately. One may also simplify the problem by considering the 

stress state in the medium as the sum of two solutions. The first is 
obtained by solving the problem for the uncracked medium under the pre- 
scribed thermal and mechanical loads, and the second is found from the 
cracked medium in which the equal and opposite of the stresses given by 
the first solution acting on the crack surfaces are the only external 

loads. In this paper the primary interest is in the stress intensity 
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factors due to transient thermal stresses. It is, therefore, sufficient 
to consider the problem under the crack surface tractions given by (19) 
and (20) (with the opposite signs). 

The plane elasticity problem shown in Fig. 1 requires the solution 


of the following equilibrium equations 

M>’ 2 ^ + 2 <^ + w> = 0 • (>1 - 2) • (21) 

9 2 Uj 3 2 V, 

( ‘r , > v2v i + 2 W + ly^’ = 0 ’ (i=1 ’ 2) • (22) 

subject to 

a lxx (0,y) - 0 , °jxy(^ ,y ^ = 0 , 0<y<«= , (23) 

a lxx (h,y) = o 2xx (h,y) , ^ (h,y) = cr 2xy (h,y) , 0<y«» , (24) 

u lx (h,y) = u 2 (h,y) , v^h.y) = v 2 (h,y) , 0<y<°* , (25) 

a lij( x ,y)- >0 » a 2i j^ x,y ^° » ( x2+ y 2 )^” » (i»j =x >y) > (26) 

<7ixy( x »°) = 0, 0<x<h; CT 2 X y( x >®) = °* h<x<0 ° » (27) 

a lyy (x,0) = p-j(x) = “O^i yy ( x , t ) , a-j<x<b-j , (28a) 

v^x.O) = 0, 0<x<a-j, b^<x<h , (28b) 

CT 2yy( x ’°) - ^2^ = ” CT 2yy^ x ’^ * ®2 <x< ^2 * (29a) 

v 2 (x,0) = 0, h<x<a 2 , b 2 <x<°°, (29b) 


where, in the usual notation u^ and v^, (k=l,2) are the x and y- 
components of the displacement vector and and a 2 -j, (i,j=x,y) are 
the stresses in materials 1 and 2, respectively. Note that y=0 is a 
plane of symmetry and, hence, the problem is considered for y>0 only. 
The stresses and displacements are related through 
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( 30 ) 


V, 3U. 3V. 

a ixx = [(K i +1) 13T + (3 " K i } IF 3 ’ (i=1 * 2) ’ 


9U . 


3V ■ 


iyy = ^^i>lF + h +1 >17 ] • (,=1 * 2) ’ 


(31) 


txy 


3U. -3V. 

= v^-gy- + -57) * Ci=1>2) . 


(32) 


The differential equations (21) and (22) are solved by expressing 
U| and v^ , (1=1,2) in terms of Fourier integrals 17], [8]. Thus, by 
defining 

<t>l ( x) = ^ v 1 ( x » + °) » * 2 ( x ) = ^< v 2 (x ’ +0) • (33) 


and by following the procedure described, for example, in [8], the prob- 
lem may be reduced to the following system of integral equations (see 
[7] for details): 

b] b 2 

| [—■ + k 11 (x,s)]<|) 1 (s)ds + j k 12 (x,s)$ 2 (s)ds 

a l a 2 

tt(1+k 1 ) 

= 4 -— ■- Pl (x) , a 1 <x<b 1 , (34) 

b l r b 2 

j k 21 (x,s)4> 1 (s)ds + | [^- + k 22 (x,s)]^ 2 (s)ds 

a i a 2 


it(1+< 2 ) 

4p 2 


P 2 (x) 


, a 2 <x<b 2 , 


subject to 




(j> 2 (x)dx = 0 


(35) 


(36a ,b) 
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if 0<a.j<b.|<a 2 <b 2 «». The kernels k - j , (i,j=l,2) are of the form 


kj j(x f s) 


oo 

* 

G. ,(x,s,a)dot , (i ,j=l ,2) 
1 J 

O 


( 37 ) 


where the functions G*. are given in Appendix A. 

* J 

To improve the accuracy in evaluating the integrals shown in (37) 
and, most importantly, to examine the singular behavior of the solution 
for cracks intersecting the interface or the boundary, (37) is rewritten 
as 


kfjCx.s) 


CO 

[G. j j(x,S,a)-G i j(x,S,a)]da 

O 



(ij=l,2) , (38) 

where G~- is the asymptotic form of G.^ for large values of a. In this 

* * J 

problem, after some lengthy but relatively straightforward analysis, it 
may be shown that (see [8] and also [7] for details) 


k ij(x,s) 

kt.(x,s) 


ktj(x,s) + k^.(x,s) , (i ,j=l ,2) , 

co 

r co 

[G ij .(x,S,a)-G.j(x,S,a)]da , (i,j = l,2) , 

o 


(39) 

(40) 


oo 

^ f oo 

j(-X*s) ~~ j G^j(XjSjOi) da ) (i *j-l »2) 9 

o 


(41) 


where the functions k.. are bounded in the closed intervals [a,,b,] and 

■ J s 

[a 2 ,b 2 ], including a-j =0 , b^h, a 2 =h and the (singular) kernels k- are 

found to be 


k 



+ k 


s2 

11 


(42) 



1 6x _ 4x 2 
s+x (s+x) 2 ~ ( s+x ) ^ 


(a 1 <(s,x)<b 1 ) , 


(43) 
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( 44 ) 


s 2 ^11 c 1? (h-x) c^q(h-x)^ 

k ll = I2h-x-s) + I2h-x-s) 2 + (2h-x-s) d 5 (a-,<(x,s)<b 1 ) , 

r = 3(m-l) _ mK 2~ K 1 = , _nM_ r = . m-1 

1 1 2(m+K.| ) 2(niK2 + 1 ) * 1 2 m+K-j * 13 m+tc-j * 

y l 

m = -1 , (45) 

c 1 d,p(x“h) 

k 12 = Fx + (s-x)^ * ( a i <x<b i* a 2 <s<b 2 ’ 


3 ( k -|+ 1 ) k - j +1 k - j +1 k - j +1 

^11 ^"('ttv+k-j ) 2(mic2 +u l ) 5 ^12 m+K«j m^+l * 

c ^21 ^22^ x ~^ 

k 21 = x^s + (x-s) z * ( a i <s<b r a 2 <x<b 2^ * 

3m(i<2 + l ) m(<2+l ) m(K2 + T) m(i<2 + l) 

^21 2(tdk2 + 1 ) + 2(iti+k-j ) * ^22 mi<£+l m+tc-j * 

,s c 2 i _ c 22 (x-h) c 23 (x-h) 2 

k 22 = (s+x-2h) + (s+x-2h) z + (s+x-2h) d s ( a 2 < ^ x » s ^ <b 2 * 


(47) 


(48) 

(49) 


(50) 


C 21 


m< 2~ lc l , 3(m-l ) r _ . m-1 r _ . m-1 ,, n 

2(m+ic-| ) 2(m<2 + T ) * 22 mi<2 + l *. 23 itik^+I 


We now observe that as long as the crack tips are away from the free 
surface and the interface (i.e., for a-|>0, b^<h and a 2 >h), the kernels 
k^. are bounded and may be treated, along with kt . , as Fredholm kernels. 
In this case the integral equations (34) and (35) have only an ordinary 
Cauchy singularity and may be solved by using any one of the standard 
techniques [9], One may also note that in the limiting cases of an edge 
crack (a^=0), a crack terminating at the interface ( b-j =h or a 2 =h), or 
the crack crossing the interface (a 2 =h=b 2 ), if the variables x and s 
approach the end points 0 or h together, then the kernels k^.(x,s) become 
unbounded and, consequently, would also contribute to the singular beha- 
vior of the solution. Such kernels are defined as generalized Cauchy 
kernels [9], 


- 9 - 



5. Singularities of the Solution 


To examine the singular behavior of the solution of integral equa- 
tions (34) and (35) we observe that (a) the left hand sides of (34) and 
(35) represent the stress components a iyy( x »°) and °2yy^ x *^ outslcle 
as well as inside the cracks in materials 1 and 2, respectively, and 
(b) the known functions p^ and p 2 , defined inside the cracks, are bounded 
but cr.jyy(x,0), (i=l,2) may be discontinuous and unbounded as x approaches 
the crack tips from outside. With these observations the singular nature 
of 4>i and <f> 2 and the stress state around the crack tips may be determined 
by a systematic application of the function theoretic method to (34) 
and (35). To do this we define the following sectionally holomorphic 
functions of the complex variable z=x+iy 


b . 


Fj(z) 




Ms) 


s-z 


ds , ( j=l ,2) 


■j 


(52) 


and express the unknown functions <f> . , (j=1,2) as [10] 

J 


<f>j(s) = 


gj (s) 


(s-ajAbj-s ) 1 * 1 


> (s j < ^ < ^j ^ * 


0<Re(y J .,e j )<l, (3=1 , 2 ) , (53) 


where the unknown functions g.(s) are Holder-continuous in the closed 

J 

intervals a j <_s<t>j , (j=l,2) and are nonzero at the end points. Following 
the asymptotic analysis given in [10], the leading terms in (52) may be 
separated and F. may be expressed as follows: 

J 


FjU) 




1 


5 s (b.) 


1 


(bj-ajlA'A (z-a/ J 


Y ^ 

(bj-aj) J sinu8 J . (z-b^ ) J 

(3=1 ,2), 


b 7 + G j (z) » 


(54) 


where the functions G.(z) are bounded everywhere except possibly at the 

J 


- 10 - 



respective end points near which they may have a weaker singularity of 
the form 


fGj(z) |< 




: jk 


(j=l,2) , 


C jl " a j * e jl < Re(Y j^ * C j2 " b j ’ E j2 < Re(e j ) * 


(55) 


In (55) k=l and k=2 correspond to end points a. and b., and A., and 

J J J K 


are positive constants. 
From the PI erne! j formula 


1 

7T 


<h(s) 


^-ds = \ CFj + (x)+Fj‘(x)] , (J-1,2) 


(56) 


and (54) it may be shown that 


b i 

J Ms) 


,v-, g.(a.)cotTry. 

A ds = J ^ i 

s-x 


g^bOcot^B. 

-' LJ J-r- + G,(x) , 


Bi Y i Y i 6 i J 

(b.-a.) J (x-a.) J (b.-a.) J (b.-x) J 


J J 


J J‘ 


a j< X< bj » (J =1 > 2 ) • ( 57 ) 


In integral equations (34) and (35) there are also other "Cauchy- 
type" integrals coming from the generalized Cauchy kernels k^(x,s). 
From ( 43 ) - ( 50) and (52) it may be seen that all these integrals can be 
expressed in terms of a Cauchy integral of the form 


1 

7 T 



*j (s) 

w o 




(j=l,2) 


(58) 


where z n is outside the corresponding branch cut and hence F_.(z) is 

^ cl J c 2 

holomorphic at z= z Q (e.g., z Q =-x<( a-|,0) in k^, z 0 =2h-x>(b-| ,h) in k^, 
z o =2h-x<(a 2 ,h) in k 22 , z Q =x<h<s in k^ 2 and z Q =x>h>s in k 21 ). From (58) 
we also observe that 
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1 

IT 


f j Us) 
Ts-Z^J 7 dS 


a j 


A F i (z o> • 7 l 


»0 


*,( S) 

TikF ds 


Td^W' 


(59) 


Since the remaining kernels kt^ are all bounded in their respective 
closed intervals, with the density functions as defined by (53) the 
corresponding integrals would also be bounded for all values of x. Thus, 
in the asymptotic analysis they may be combined with the other bounded 
terms p-j and p 2 . 

The singular behavior of <f>i and <f> 2 may now be determined for any 

location of the crack tips (i.e., for a-|>0, b-|j<h, a 2 >_h) by substituting 

from ( 53)- (59) into (34) and (35), and multiplying both sides of the 

equations by (x-a.) J , or (b.-x) J, and letting x-*a., x->b. for a,>0, 

J Y1 Jgi yp J J 1 

b^<h, a 2 >h, and by x , (h-x) , (x-h) and letting x->0, x-*h-0, x=h+0 

for the limiting crack tip locations a-|=0, b-j =h , a 2 =h. In the standard 

case of cracks embedded into homogeneous materials, that is for a,>0, 

c *p 

b,<h, a 0 >h, k. . as well as k.. are bounded and the only contribution to 

I ^ * J ^ J ^ 1 

the singularity comes from the Cauchy kernels (s-x) in (34) and (35). 
Thus, since g-(a.)?*0, g-(b.)/0, from (57), (55), (53), (34) and (35) it 

J J J J 

may easily be seen that 

cotiTYj = 0, cotTrgj = 0; Yj = 1/2, 6j = 1/2 , (60) 

leading to the standard results for the embedded cracks, namely 


g.-(x) 

♦ ,(x) = * , (a .<x<b. ) , (j=l ,2) . (61) 

The characteristic equations for y,- and e- corresponding to various 

J J 

other special crack geometries may be obtained by using the asymptotic 
expressions given in this section and the related generalized Cauchy 
kernels k*j(x,s). If, for example a^O (i.e., the case of an edge crack), 
then following the procedure outlined above, substituting from (43), 

(54), (57) (with a^=0) and (59) into (34), multiplying both sides by 
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Y-i 

x , letting x-* 0 , and by observing that g^O^O, the characteristic 
equation for y-j is obtained to be 


COStty-j 


2(y r l) 2 + 1=0. 


( 62 ) 


The only acceptable root of (62) is y-|= 0, that is, as expected at x=0 
the derivative of the crack surface displacement is found to be bounded. 
If b^<h, a 2 >h, then the results found for , y 2 and 62 in (60) would, 
of course, remain valid. 

Consider now the case of a crack tip terminating at the interface. 
Let, for example, b^=h, a^>0, a 2 >h. Clearly for the embedded crack 
(a 2 ,b 2 ) (35) would again give y 2 =b 2 =-1/2. Also, regardless of the loca- 
tion of the crack tip x=b-j , (34) would give y^=-l/2 for a-|>0 and y^=0 

for a,=0. Thus, to determine e, (defined in (53) for b-,=h) it is suffi- 

-1 s2 

cient to consider only the singular kernels (s-x) and k^(x,s) in 
(34). By substituting from ( 57)-( 59) (with b^=h, z 0 =2h-x) into (34) and 
by separating the terms that are singular at x=b^=h, it may be shown 
that 


g 1 (h)cot^ 6 1 1 g-j (h) 

- — - — - 

(h-a.|) (h-x) (h-a-j) sinirg^ 


[ c n +c 12 (h-x) 


_d_ 

dx 


+ j c 13 (h-x ) 2 ^-] — y- = H(x) , a-j<x<h , (63) 

(h-x) 1 

where H(x) represents all other terms in (34) that are bounded at x=h. 
After carrying out the differentiation, if we multiply both sides in 
(63) by (h-x) 1, let x+h and observe that g-|(h)^0, we find 

costtS-j “7) c i3&i($i + l ) ” c i 2^1 ” ^11 ~ ^ » (64) 

where the coefficients c-jj are given by (45). 

Similarly, for relative crack tip locations 0<a^<b^<h = a 2 <b 2 we 

find 
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(65) 


*J 

costty 2 + j c 23 y 2^ y 2 +1 ^ + c 22 y 2 + C 21 = 0 

with c 2j . as given by (51). 

Finally, to examine the singular behavior of and <j> 2 at x=h for 
a crack crossing the interface we let 0<a^<b 1 =h=a 2 <b 2 and observe that 
the end point x=h is common to both cracks. In (53) we must, therefore, 
have 


B 1 = Y 2 = 0<Re(e)<l, (b^=h=a 2 ) . (66) 

In this case by separating the terms that are bounded at the common end 
point x=h the integral equations (34) and (35) may be expressed as 


1 

TT 


h 

C^-+ k^(x,s)]<t>-| (s)ds +1 


b ? 

( c. 

k 12 (x,s)<f> 2 (s)ds = H^x) 
h 

a-j<x<h , 


9 


1 

7T 


h 


k 21 (x,s)ct>.,(s)ds +1 


[g_^ k 22 (x,s)]c}> 2 (s)ds = H 2 (x) , 

h 

h<x<b 2 , 


(67) 


( 68 ) 


where H-j and H 2 represent all other terms in (34) and (35) bounded at 
x=h. By substituting now from (57) with b-|=h, a 2 =h, B-|=y 2 = B> and (54), 
(58) and (59) with z Q =2h-x and z Q =x into (67) and (68) we obtain 


g^h) 


-g,(h)cotiTB 

— 1 + 

Y i B Y 1 

( h-a-j ) (h-x) ( h-a-j ) simrB 


[ c n+ c l 2 (h- x ) d d x + 2 3 (h-x) 


g 2 ( h ) 


Bo L 12 
(b 2 ~h) sinuB 


[d 1 o(x-h) ^ + d 1 - - = H-j (x) , a ] <x<h, (69) 


(h-x) 


B 


-9-1 (h) 

Y i 

(h-a-j) sinirB 
„ 9 2 (h) 


[d 22 (x-h) Hv - d 91 ] 


dx 2V 


U-h) 1 


g 2 (h)cot7rB 

Bp 

(bo-h) 2 


(x-h) f 


(bp-h) sinTTg 


[c 21 -c 22 (x-h)i+^ < x - h ) 2 ^7 J rr6 = H 2 (x) ' h<x<b 2' 


dx 


dx - 1 


(x-h) 


B 


(70) 
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Multiplying (69) and (70) respectively by (h-x) e and (x-h) 6 and letting 
x-*h, it may be shown that 


i 9i(h) 

[c-| i 2 2 ^ ( 1 6 ) B~costt 3] 


Y 1 

( h-a-j ) simrg 


+ (^11 ~^i 2^ ) 


g 2 (h) 


= 0 


2 . 


( b^-h ) sin-rrB 


( d 2 ^j +^ 22 ^) 


g-, (h) 


Y*i 

(h-a-|) sinuB 


1 g 2 ( h ) 

+ [coS7rB+C2i + BC22 + 2 3 C 3 + l ) c 23-^ g 0 • (71 a, b) 

(b2"h) ^simrB 

Since g-j(h) and g2(h) are nonzero, from (71) we obtain the following 
characteristic equation to determine 8 and the relationship between g-j(h) 
and g 2 (h): 

C-.0 c o3 

[c-| 1 +C 1 2® + ~~t~ 6 n + B)-COS 7 T 6 ][c 21 +C 2 2 8 +— 8(l+6)+COSirB] 

-(d 11 -d 1 28)(d 2l +d228) = 0 , 0<Re(8)<l , (72) 

Yi 1 

9-j(h) (h-a^) c 2i +c 22 e + 2 c 23 8(1+8 )+costtB 

g 2 ^ (b _ h) 3 2 d 21 +d 22 e * (?3) 

A close examination of the characteristic equations (64), (65) and (72) 
would show that in each case the acceptable leading root is real for all 
material combinations. 

It may be remarked that for a-|>0 the solution of the integral equa- 
tions (34) and (35) contain two arbitrary constants [10], [9]. In the 
case of two non-intersecting embedded cracks (that may include b-|=h or 
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a 2 =h) the single-valuedness conditions (36) are used to determine these 
constants. For a^=0, only (35) contains an arbitrary constant 

which is determined from (36b). For the crack crossing the interface 
(36) reduce to the following single condition which accounts for only 
one of the constants: 

h b 2 

| <f>-j(x)dx + j <j> 2 (x)dx = 0 . (74) 

a^ h 

The additional condition needed to determine the second constant is pro- 
vided by (73). 


6. Stress Intensity Factors 


Recalling that the left hand sides of equations (34) and (35) rep- 
resent aiyy(x,0) and °2yy^ x ’^ outside as well as inside the cracks, 
from (34) and (35) we can write 


*| Hhf<£ 

"VT a iyy (x,0) 


1 

7T 


. b s 




+ k! j .(x,s)] 4 > j .(s)ds + M.(x), (i=l,2) 


a . 
J 


(75a, b) 


where functions and M 2 correspond to the terms that are bounded at the 
end points. If the cracks are embedded into the respective homogeneous 
materials (i.e., if 0<a-j<b^<h<a 2 <b 2 ), then the kernels k^. are bounded 
for all values of x and s making also only bounded contributions to 
a iyy( X *0) . * n tbls case we °b serve that for x 0L, tside the cracks the 
Cauchy integrals in (75) may be expressed by (58) and (54), which, upon 
separating the singular terms, and substituting into (74) give 


4p k 9t c ( a b) i 
o k yy( x »°) " 1+ I 


9k (b k> 


k /b k -a k /x-a k /b k ~a k /x-b k 


+ P k (x)], (k=l ,2) (76) 


where P-j and P 2 represent all the bounded terms on the left hand side 
of (34) and (35). Defining now the mode I stress intensity factors by 
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(77) 


k,(a.) = lim 
1 3 x->-a j -0 


/2(aj-x) a j yy ( x >°)> (J =1 » 2 ) 


k,(b.) = lim a ivv (x,0) ,(>1,2) 

J x-»-bj+0 J Jyy 

from (75) it is seen that 


(78) 


4v i g,(aj 

1/ ( a \ s J J J 

r > 1 +k, 


j *9977 


• MV ■ 


Uk s 7^997? 


(79) 


We now consider the crack terminating at the interface and let 
°l a l» b -| =h » h<a 2 <t> 2 . Again, let the stress intensity factor at x=h be 
defined in terms of the local cleavage stress, which in this case is 
a 2 yy(x,0), x>h. Thus, analogous to (78) we may define 


k,(h) = lim n (x-h) cr ?hh (x,0) . 
1 x+h+0 


(80) 


On the other hand, keeping only the terms that may contribute to the 
stress singularity and by substituting from (48), (75b) may be expressed 
as 


l+< 2 

4y 2 


a 2 yy(X,0) l 



d 9 o(x-h) 

) 4>i ( s)ds + Q-j(x) , 


(81) 


where Q-|(x) represent all the remaining terms that are bounded at x=h. 
Since x is outside the cut (a^,h), (58), (59) and (54) may be used in 
(81) with z Q =x, giving 


l+« 2 

a 2yy (x>0) 


-d 2 -|F^(x) + d 22 (x-h) F-j (x) + Q-|(x) 


(82) 


or 

a 2yy( x> °) 


4v 2 g-j (h) 

1 + ko Ti 

(h-a-j) siriTte-j 


( d 21 + d 22^1 ^ 


(x-h) 6 ' 


+ Q 2 (x) 


(83) 
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f Jc\ 

where Q again is bounded or, at most may have a singularity weaker v ' 
than (x-h ) . From (83) it is seen that the cleavage stress c^yy has 

a singularity at x=h with power b-| and from (80) and (83) it follows 
that 

4y ? SI g-, (h) 

¥»> ■ vS, < d 21 +d 22 S l> < 84 > 

c ( h-a-j ) sinu^i 

where B-| is the acceptable leading root of (64), g^ is defined by (53), 
d 21 and d 22 are given by (49), and y-|=1/2 for a.j>0 and y^=0 for a^=0. 

Similarly, for the crack in material 2 and terminating at the inter- 
face, i.e., for 0<a.|<b.|<h, h=a 2 <b 2 , from the asymptotic examination of 
<7lyy(x,0), x<h we find 


k,(h) = lim SI (h-x) Y2 o, w (x,0) 
1 x->h-0 lyy 


4y 


1 


SI g 2 (h) 


1+K 


1 P? 

(b 2 ~h) sinTTY 2 


(d ir d i2’'2 ) 


(85) 


where y 2 1S given by (65), d^ and d^ 2 by (47) and b 2 =1/2. 

In the case of a crack crossing the interface the stresses as well 

as the displacement derivatives would have a singularity at x=h of the 

form r“ e , r being the radial distance from the point x=h, y=0 and b the 

root of (72). In this problem of particular interest are the normal 

and shear cleavage stresses a xx and along the bond line x=h. The 

expressions for cr VY (h,y) and a Y ..(h,y) may be obtained in terms of 
xx xy (**) 

4>1 and <f> 2 by going back to the main formulation of the problem' 

J*1 This is due to the possible singular terms coming from G-j(z) that 

appear in (54) and has the behavior as given by (55), and from the 
fact that the characteristic equation (64) may have more than one 
root satisfying 0<Re(B-j)<l. 

(**) Note that, because of continuity we have a -| xx ( h >y) = ^2xx^ h,y ^ =a xx^ h,y ^ 
and a ixv^ h5y ^ G 2xy (h ’ y ^ =a xy^ h,y ^ 
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After somewhat lengthy analysis and upon separating the singular parts 
of the kernels it can be shown that (see [7] for details) 


o xx (h,y) 


°xy (h,y) 


Zv Z l 1 

l+<2 1 IT 


Z 2z 1 1 

1 +k 2 i ^ 


C k 3j( y *s) + k 3j Cy»s)]4>j(s)ds, (0<y) , 

a j 

b j 

j (y»s) + k£.(y,s)]<j> J .(s)ds, (0<y) , 


( 86 ) 


(87) 


f 

where b-^a^h, the kernels k..., (i=3,4; j=l,2) are of the Fredholm type 
and the generalized Cauchy kernels k^., (i=3,4; j=l,2) are given in Appen- 
dix B. In (86) and (87) only kL contribute to the stress singularity 

* J 

at (y=0, x=h). The singular terms in (86) and (87) may be evaluated by 
using again (58), (59) and (54). For example, by substituting from 


g-. (s) 

^(s) = £ , (a 1 <s<h) 

(s-a-j ) ] (h-s) 


( 88 ) 


into (86), using the equations (Bl), (54) and (58), the first term in 
(86) may be obtained as 


1 

7T 


4>1 (s)ds 

(h-s'FV 


1 1 

2 TT 


T I 

^h-s+iy + h-s-iy^l ( s ^ ds 


'1 


g n ( h) , 

T +R(y) 


2( h-a-j ) Y] sin ^ y 


(89) 


where R represents the terms that are bounded at (y=0, x=h). It is then 
clear that the leading terms in the asymptotic expressions for a (h,y) 

AA 

and a (h,y) will be of the form y , where g is the root of (72). We 
xy 

may now define the following stress intensity factors to characterize the 
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stress singularity at (y=0, x=h): 


k xx = 1<fm yBa xx^ h » y ^» k xv = lim y$a xv^ h » y ) • (90a, b) 

xx y ^+0 xx xy y ^+0 xy 

By using the results of the form (89), from (86), (87) and (88) it may 
then be shown that 


k 


xx 


"2 

(1+K 2 )sin ^ 


j- A+(l-2g)B 

(h- ai > Vl 


»,(h) _B±I3=2filA , 

(b 2 -a) 2 


(91) 


W (i +K2 )cos f 


( ft ~ ( - 1 ~ 26 ! E g,(h) + 6 ~ (1 ' 2B > A g,(h)] 


(h-a-j ) 


y l 


(92) 


(b 9 -h) 


A = 


m( 1 +^2 ) 
mK 2 +l 


m(< 2 +l ) 

B - _ . » m - u i /y 2 


m+tc 


1 


(93) 


where e 2 =l/2, y-j =1 /2 for a^>0 and y-|=0 for a-j=0. 

Two stress intensity factors defined by (90) at the same point is 
simply for convenience. As one might expect and as seen from (73), (91) 
and (92), g-|(h) and g 2 (h) are not independent and there is only one 
parameter characterizing the stress singularity. If (r,e) are the polar 
coordinates with the origin at the point (x=h, y=0), for small values 
of r the stresses may be expressed as 

k 

o^-fr.e) ■-fff(e) , (i »j=x,y) , 0<e<n- (94) 

«J pP IJ 

Thus, from (90) and (94) it follows that 

k xx ■ k o f xxW 2 > • k xy ■ k o f xy'" /2 > ’ < 95 > 
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7. On the Solution of the Integral Equations 

For numerical solution the integral equations (34) and (35) are 
first normalized to read 

1 6 . . 1+K . 

7 j C—J- + m f j( r »p)]^j(p)dp = -^ J -q i (r) , -l<r<l, (i=l,2) (96) 

-1 


where 


r = 


b_.+a. „ b..+a.. 


»=I^(S-V) . q,(r) ■ Pi(x). 


1 1 


1 1 


b .-a . 

m ij( r,p ) * \ J k ij( x » s )» ♦ j (pH J (s)- 


Vp) 


Y-; e,-* 

(i+p) j O-p) 3 


(97) 


The unknown functions f-j and f 2 are real and bounded and may be determined 
by using any one of the techniques available for solving singular integral 
equations [9]. In this paper it is assumed that 

N j 

fj(p) = S A jn p n , - 1< P <1 , (j=l , 2 ) . ( 98 ) 

The unknown coefficients are then determined by using the special formu- 
las and the procedure developed in 111] and [12] for various crack con- 
figurestions. 


8. Results and Discussion 

As examples for the thermal shock problem described in the previous 
sections we consider two material combinations. Material pair A cor- 
responds to a stainless steel layer welded on a much thicker ferritic 
steel medium approximating cladded pressure vessels of very large dia- 
meters. Material pair B represents a ceramic layer (Material 1) bonded 
to a thick steel substrate (Material 2). The properties of the two 
material combinations are given in Table 1. Since the problem is for- 
mulated in terms of dimensionless quantities, it is sufficient to con- 
sider only the ratios shown in the table. 
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Table 1. Properties of material pairs used in examples 


• < 



¥ k i 

Dg/D] 

o^/oii 

e 2 /e, 

v 2 /v l 

Material pair A 

3 

3 

0.75 

1 

1 

Material pair B 

3.385 

4.070 

2.2939 

0.6111 

1 


Figures 2 and 3 show some sample results for the normalized stresses 

a iyy/ 0 o defined by (19) and (20) in the material pair A. Note that these 

are essentially the plots of normalized temperatures (see equations (16)- 

(20)). Comparison of the two figures clearly shows the influence of 

ramp cooling on the transient thermal stresses, particularly for small 
(*) 

values of time v ' . Stress intensity factors for various crack geometries 
in the material pair A are shown in Figures 4-16 (see Fig. 1 for nota- 
tion). In Figures 4-16 all stress intensity factors are normalized with 
respect to 


k o = V* ’ °o = 


: l a l ( 


(99) 


where l is the total crack length for an edge crack (i.e., £=b-| for a-j=0, 
b^<h, and l=b ^ for a-|=0, b-ph 5 ^) and the half crack length for an inter- 
nal crack. The parameter t q defined by 


r 0 = yyh 2 (ioo) 

is the measure of the rate of change of surface temperature during ther- 
mal shock, t being the actual duration of the ramp (see (6a) and Fig. 
lb). 

Figure 4 shows the variation of the normalized stress intensity 
factor k 1 (b 1 )/cr 0 v^ 1 for an edge crack and for t o =0 with the Fourier 


v 1 In the materials used for a clad thickness h=6 mm. t=10r, t and x 
being the real and normalized times (in units of seconds) (see 
(18)). 
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number r=tD-|/h 2 representing the normalized time. Note that in the 
material pair A E-j =E2 and v-|=v 2 . Consequently as t-*° k^/k 0 would 

approach the uniformly loaded half plane solution of 1.1215 (see, for 
example [11]). This asymptotic trend may joe seen in Fig. 4. For a 
fixed crack length b.j=0.5h the influence of the cooling rate on the 
surface x=0 as measured by x Q is shown in Fig. 5. It is seen that for 
small values of time this influence could be quite considerable. 

The stress intensity factor given in Fig. 4 as well as that given 
in all subsequent figures are presented in normalized form and the nor- 
malizing stress intensity, k Q always contains a length parameter (as, 
for example, in k 0 =CT 0 vf>^ for the case shown in Fig. 4, b-j being the crack 
length). Thus, in considering questions relating to fracture propaga- 
tion and fracture instability, the normalized stress intensities must 
be multiplied by k Q before they are compared with the fracture tough- 
ness, Kj£ or used in a subcritical crack growth model. One consequence 
of this in connection with the example described in Fig. 4 is that 
k^ ( b^ ) is a monotonously increasing function of both crack length and 
time implying an unstable fracture process. 

For an edge crack crossing the interface the results similar to 
Figures 4 and 5 are shown in Figures 6 and 7. For each crack length 
the figures also show the asymptotic values of the stress intensity 
factor as t-*». The stress intensity factors k^(a 2 ) and k-|(b 2 ) for an 
under-clad crack are shown in Figures 8-11. Again, the figures also 
show the asymptotic values of k-j as t-*=°. One may note that in material 
pair A, despite the large differences in thermal coefficients a, D and 
k', the elastic constants E and v for the two materials are the same. 
Consequently, for a crack terminating at the interface we have y 2 =1/2 
(Figures 8-11), g-|=l/2 (Fig. 12) and for a crack crossing the interface 
8=0 (Figures 6,7,13-16). Figure 12 shows the stress intensity factor 
for an edge crack terminating at the interface (b-|=h). Figures 13-16 
show the results for material pair A containing an internal crack cross- 
ing the interface (0<ai<b^=h=a 2 <b 2 «»). Figures 13 and 14 show the 
influence of ramp duration r o on the stress intensity factors for a 
fixed crack geometry, a^=0.2h, b 2 =2h. Figures 15 and 16 gives the stress 
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intensity factors for an under-clad crack growing into the clad. The 
figures also show the asymptotic values of the stress intensity factors 
for t-*®. 

The results for the material pair B are given in Figures 17-25. 

Figures 17 and 18 show some sample results for the transient thermal 
T T 

stresses tf-jyyCx**) and ° r 2yy^ x,t ^ obtainecl ‘f or T o =0 and T o =60 * Again, 
from (16) and (17) it follows that these are essentially the temperature 
distributions in the composite medium. For these materials the stress 
discontinuity at the interface is given by (see (19) and (20)) 


E 2 a 2^’ v l^ 

crjy/h.t) W , " v 2 ) 


1.4019 


(101) 


Figures 19 and 20 show the normalized stress intensity factor for 
an edge crack in material 1. The limiting case for the edge crack ter- 
minating at the interface is shown in Fig. 21. In this case for the 
material combination under consideration (64) gives the power of singu- 
larity as e-j =0.552538 . The normalizing stress intensity factor used in 
Fig. 21 is defined by 


k o = ^ h S] (102) 

The results for an edge crack going through the interface (i.e., 
a^=0, bi=h=a 2 <b 2 <« , J Fig. 1) are shown in Figures 22-25. Figure 22 shows 
the effect of the ramp duration t on the crack tip stress intensity 
factor ki ( ) for a fixed crack length b2=2h. The influence of the 
crack length b 2 on k^(b 2 ) is shown in Fig. 23. Finally, Figures 24 and 
25 show the stress intensity factors k and k governing tensile and 

XX AJ 

shear cleavage stresses at interface, (x=h, y>0). 
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APPENDIX A 



H 3 = d 2 d 3 [l+2ah(2at-1)]e" a(t+2h) -d 1 d3e' a(t+4h) 

-d 2 d 3 [1+2a(t-h)]e a(t “ 2h ^+d 1 d 3 (l+2at)e a(t ' 4h ^ , 
H 4 = d 2 d 5 e~ 0£t -{d 1 d 5 +2d 3 d 5 ah[l-2 a (t-h)]}e‘“ (t+2h) , 


H 5 2^ e “ D^ H 3 “ k 1 H 7^ * 

H 6 = ^ (k 1 H 8“ H 4 ) 1 

H 7 = d 2 d 4 (l-2at)e" at - d 2 d 3 [l-2 a (t-h)]e' a(t+2h) 

+ {2d 2 d 3 ah[l+2 a (t-h)-d 1 d 4 }e a(t " 2a, +d 1 d 3 e a(t " 4h) , 
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H 8 = {d 4 d 5 [T“2a(t-h)-2d 2 d 5 ah}e" at 


- d 3 d 5 n-2a(t-h)]e" a(t+2h) , 


H 9 = ]t(K 2 -2a1)H n -^H 3 e 2 » h ] 

H 10 " 'jC('=2" 2 “ h ' H 12 " Dijj H 4 e 


» 


do 


* 


H n ■ £ < e0t + n V • 


D = -d 2 d 4 +[d 1 d 4 +d 2 d 3 n+4h 2 a 2 )]e‘ 2ah -d 1 d 3 e‘ 4ah , 

d-| ~ niK^-Ki 9 d 2 - m< 2 +l 9 d 3 - m— 1 9 d 4 - 9 dg 

dg = m(vci+1), m = y-j/p 2 • 


(A5-A16) 

(A17) 

(A18) 
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APPENDIX B 


The generalized Cauchy kernels that appear in equations (86) and (87). 



(Bl) 

k 32 (y,s) - • 

(B2) 


(B3) 

fc s _ By + Ay[y 2 -3(s-h) 2 ] 

k 42 (y,s) ( s-y ) z +y z L(s-y)*+y*J z * 

(B4) 

m( <2 + 1 ) mUg+l ) 

(B5) 

A . l 9 B , 9 IT] p -i l U r\ 

m<2 + l m+K-j 1 2 
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APPENDIX C 
List of Symbols 

Crack tip locations, i=l,2 
Specific heat, i=l ,2 

Coefficient of thermal diffusivity, i=l,2 
Modulus of elasticity, i=l,2 
Thickness of the layer 

Stress intensity factors at the crack tips, i=l,2 
Normalizing stress intensity factor 
Coefficient of heat conduction, i=l,2 
time 

Ramp duration 

Surface temperature for t>t Q 
Initial temperature 
Temperature, i=l,2 

Components of the displacement vector, i=l,2 
Rectangular coordinates 
Coefficient of thermal expansion, i=l,2 
Poisson's ratio, i=l,2 

Elastic constants, y=E/2(l+v), k=3-4v (plane strain) 
Mass density 

Normalizing stress, a o =-Eiai 0 o /(l-vi ) 

Stress components, (i,j)=(x,y) 

Temperature, e-'T.-T*,, i=l,2; e^^-T* 

Normalized time (Fourier Number), T=D-|t/h 2 
Ramp duration, T 0 =D.jt o /h 2 


y 



(a) 


Ftg, 1 The crack geometry and the ten 




.20 .40 .60 .80 1.00 


f 



material pair A) 


Fig. 2 



Fig. 3 






in 

cu 



r-tDJh* 


Fig. 5 The influence of the ramp duration t q on the normalized stress 
intensity factor k-|(b 1 )/ 0 o /T 1 for a fixed crack length b-j/h= 
0.5; U-pb-j), a.j=0, material pair A) 
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Fig, 6 Normalized stress intensity factor kifbgj/o /£ in the medium 
with an edge crack crossing the interface for t q = 0; (£=b 2 > 
a-j=0, b^=h=a2, material pair A) 
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Fig. 7 The influence of surface cooling rate or t q on the normalized 
stress intensity factor k-|(b 2 )fa Q SZ for an edge crack cross- 
ing the interface; (f=b 2 , b 2 /h=4, material pair A) 
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Fig. 12 The normalized stress intensity factor k(b.|)/a 0 /Ej" for 
broken clad; (a^O, £-|=h, f>2 =a 2* matena ^ pair A) 












O 

O 


O 

cu 

cn 


o 

^r 

cu 


o 

CO 


o 

CO 


o 

o 


r — * oo 

'/* *(«,)/ T/v/i/J 



.0 B.O 16.0 24.0 32.0X10 

T = tD x ll i 2 


Fig. 15 The normalized stress intensity factor k-|(a ^)/o Q /UT for an 
embedded crack crossing the interface; C^ s= E> 2 “ a l * 
t o = 0, material pair A) 
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Ffg. 16 The normalized stress intensity factor k 1 (h 2 )/a 0 v^7^ for an 
embedded crack crossing the interface; U=b 2 -a-|, b 2 =2h, 
t o =0, material pair A) 
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Fig. 20 The influence of t 0 on the normalized stress intensity factor 
k-|Cb-|)/ 0 o v^j" in material pair B with an edge crack; (£j s b^» 
a-|=0, b-j/h^O.B, 
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>n the normalized stress intensity factor 
rial pair B with a broken layer; (a-j=0, 
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Fig. 22 The influence of t q on the normalized stress intensity factor 
^•|(b 2 )/c 0 /E in material pair B with an edge crack crossing 
the interface; (a^=0, b^h^* b 2 = 2 h, o 0 =-E^a-j6 0 /Cl-v-j)) 
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The normalized stress Intensity factor in material 

pair B with an edge crack crossing the interface; (t o =0, 
l=b 2 , a,=0, 5.=h=a ? , a. r r"*, 6=0.01872238) 
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The normalized stress intensity factor in material 

pair B with an edge crack crossing the interface; (t q s 0, 
l=b 2 , a^=0, b 1 =h=a 2 , 6=0.01872238) 


